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1 Introduction 

The nature of dark matter (DM) is still a mystery in our knowledge about the universe. In 
spite of increasingly accurate cosmological and astrophysical observations confirming the 
existence of cold DM (see e.g. [1, 2] for recent reviews), its properties, such as spin, mass 
and flavor structure, remain elusive. The most often considered theoretical candidate for 
cold DM is a weakly-inter acting massive particle (WIMP), which may be detected via non- 
gravitational interactions. However, there are no candidates for WIMPs in the standard 
model (SM) of particle physics. Therefore, any theory trying to describe DM requires an 
extension of the SM. 

Searches for DM have been carried out in many ways, including cosmology and astro¬ 
physics observations, direct- and indirect-detection experiments, and searches at particle 
colliders. Recently there has been a rise in the interest in collider searches for DM, as 
the LHC opens up new energy regimes and more search possibilities, and its sensitivity 
continues to be improved as the integrated luminosity increases. The generic DM signal at 
the LHC consists of the production of one or more SM particles accompanied by missing 
energy. Using the data accumulated in the first stage of running, the ATLAS and CMS 
Collaborations have already searched for DM in events with missing transverse momentum 
associated with a single particle X, i.e. a mono-A signal, where X may be a jet, a photon. 
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or a W or Z boson [3-8]. So far no excess has been discovered, and limits have been set 
on the energy scale of the DM interaction and the mass of the DM particles. 

Stringent constraints on spin-independent and spin-dependent DM interactions have 
already been obtained, respectively, by the LUX experiment [9] and by mono-X searches 
at the LHC [3-8] . These limits have been derived in the framework of effective field theory 
(EFT) [10-14], which has the advantage of being model independent and involving only 
a few parameters, thus allowing for easy comparison between theoretical predictions and 
various kinds of experiments. However, the EFT is only a low-energy approximation of 
more complete models. Its validity should be carefully checked when applied to LHC 
processes. The issue has been widely discussed in the literature, see e.g. [11, 12, 14-19]. A 
potential problem is that the requirement that the fundamental energy scale of the EFT 
be much larger than the momentum transferred to the DM particles is not always satisfied 
over the entire region of phase space probed by the experimental searches at the LHC. For 
example, a mediator with mass in the few-TeV range, which one would integrate out in the 
EFT approach, could become accessible in events with very high px, thus rendering the 
EFT assumption invalid. This issue will become more pressing in Run-H of the LHC. As 
a result, an alternative framework based on so-called simplified models has been proposed 
as the standard for future DM searches [20-38] (see [39, 40] for recent reviews). The 
idea behind simplified models is to incorporate the most relevant degrees of freedom of 
the underlying (UV complete) model at accessible energies, while keeping much of the 
simplicity of EFT. Moreover, they may provide new search signals besides the missing- 
energy signature, such as the direct production of the mediators. Meanwhile, simplified 
models also add some complexity. Unlike in EFT, the couplings between the mediator and 
the DM and SM particles need to be specified, and the width of the mediator, which enters 
as a new parameter, needs to be calculated. As a consequence, the results can depend on 
the couplings in a rather nontrivial way. 

Current theoretical predictions for DM production at the LHC have mostly been pro¬ 
vided at leading order (LO) in QCD, although some next-to-leading order (NLO) QCD 
results exist within the EFT framework [41-46]. With the data of the 13 TeV collision, 
it is expected that more precise constraints can be obtained. In order to match the ex¬ 
perimental accuracy, it is important to have complete theoretical predictions with NLO 
accuracy, including parton-shower effects. Such higher-order corrections tend to increase 
the cross sections by a few tens of percent and reduce the scale uncertainties, thus mak¬ 
ing the theoretical predictions more reliable. In addition, with more final state particles, 
kinematic distributions can be predicted more reliably at NLO, serving as a reference for 
experimentalists to set appropriate kinematic cuts. 

In this paper we provide theoretical predictions for mono-Z production at the LHC 
including NLO QCD corrections and the parton shower. To the best of our knowledge, 
such predictions have not been presented before. The mono-Z channel is particularly 
important for DM searches at the LHC. Previous studies of this channel within the EFT 
framework have been based on LO predictions [47-51]. The ATLAS Collaboration has 
reported constraints on the pp ^ Z -|-^t cross section, where the Z boson is reconstructed 
either via a leptonic decay [52] or a hadronic decay [5]. Our results provide a more solid 
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Figure 1. Illustration of the connection between SM and DM particles via messengers Yi. 

theoretical basis for future studies in this channel. 

We emphasize that our work is part of a larger project aiming at providing complete 
NLO predictions for all relevant DM search channels at the LHC based on the automated 
framework of FeynRules [53] and MadGraph5_aMC@NLO [54] . As a first step, we have 
implemented the simplified models with spin -0 and spin -1 s-channel mediators, shown in 
figure 1. This work has been coordinated with the authors of [55]; however, compared 
to this reference we add some specific interactions and corresponding higher-dimensional 
operators contributing to the V -j- production channels (with V = Z,W, 7 ) so as to 
allow for more complete studies of mono-H signatures. There are two main features of 
this approach. First, it realizes full automatization from drawing Feynman diagrams to 
generating collision events at NLO matched with parton shower and is thus a convenient 
tool for both theorists and experimentalists. Second, with a few extra lines of code one 
can in principle calculate any mono-W signal, and hence complete NLO results for DM 
production become readily available. In this work, we focus on the mono-Z production to 
illustrate the capabilities of such an approach. Throughout this paper we stick closely to 
the benchmark points suggested by the ATLAS/CMS DM forum [40], but we stress that 
our implementation is far more general and different scenarios can be studied as well. An 
analogous application of this framework to the mono-jet channel has recently been discussed 
in [55]. In addition, some loop-induced DM production modes have been explored in [56] 
based on recent tools for calculating loop-level processes [57] . The analysis of other search 
channels is left for future work. 

This paper is organized as follows. We first introduce the simplified DM models and 
effective operators in section 2 and describe their implementation in section 3. We then 
present QCD NLO results for total cross sections in section 4 and kinematic distributions 
in section 5. The discovery potential of mono-Z signals at the 13 TeV LHC is discussed in 
section 6 . Section 7 contains some conclusions and outlook. 

2 Simplified models for dark matter production via s-channel mediators 

We start by defining the simplified models underlying our analysis. We consider three types 
of DM particles W, namely a Dirac fermion Xd, a real scalar Xji, and a complex scalar 
Xc- We also consider two kinds of mediators Tj, namely a real scalar Yq and a vector li. 
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see figure 1. These particles are assumed to be singlets under the SM gauge group. Our 
notations for particles and coupling constants are consistent with those employed in [55], 
but we explore a larger number of interactions in the present work. Model files, which 
specify our simplified models, have been released on the website [58]. 

We consider different types of interactions between the SM, the dark sector and the 
mediators. The renormalizable interactions of the messengers with the DM are described 
by the Lagrangians 


= o [Xhid.Xc) - {d^X*c)Xc]Yl^ + XDl,{gl, + 9L^^)Xd , 

1 _ ( 2 . 1 ) 

4 m = ^Mx^glc^XRXRYo + Mxag^XcXhXcYo + XD{g^Xo+i9^o^5)XDYo. 

We omit quadratic terms in the messenger fields, which will play no role in our analysis. 
The leading interactions of the scalar messenger Iq with scalar DM fields are described by 
dimension-3 operators. We use the DM masses to serve as natural scales in these terms, 
but one is always free to readjust the coupling strengths if a different normalization scale 
is employed. The renormalizable interactions of the messenger fields with SM quarks are 
described by the Lagrangians 


+ 9i^l5)dj + uaf,{g^^. + gu,^75)uj 


-'1 > 




-V 


di % ( 4 , + ^9L75)dj +Ui^ {gl. + ^ 4 , 75 ) 


rj.U 

Vi /_s 


v/ 2 ' 


^/ 2 ' 


Ui 


( 2 . 2 ) 


Tn 


which refer to the mass basis of the SM quarks. In order to avoid large flavor-changing 
neutral current interactions, we will assume that the coefficients g^.. are flavor diagonal, 
i.e. = dij 4 B = V, A, S, P. We also assume that the (pseudo-)scalar couplings are 
aligned with the SM Yukawa interactions, and hence in the mass basis they are proportional 
to the Yukawa couplings yf of the various quarks. It then follows that the interactions of 
the scalar mediator Yq with light quarks are strongly suppressed. As a consequence, this 
type of coupling is best probed in the search channel tt+^T [55, 59, 60]. We emphasize that 
the operators containing Yq in (2.2) are gauge invariant only after electroweak symmetry 
breaking, which implies that the coefficients gq^. must implicitly contain a factor u/A [v is 
the vacuum expectation value of the Higgs field and A is a new physics scale), suppressing 
these interactions further even for the case of the top quark. 

The Lagrangians (2.1) and (2.2) have been discussed previously in [55]; we show them 
here only for completeness. In addition, we also consider couplings of the mediators with 
the electroweak bosons Z, W and h of the SM. For the Yi mediator, the corresponding 
renormalizable couplings are derived from the Lagrangian 

4w = 9^1 - (D^<?i)t</.] Yf , (2.3) 

where (p is the Higgs doublet. A kinetic-mixing term with the hypercharge field of the form 
d^Yp' can be reduced to the operators contained in in (2.2) and in (2.3) 
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using an integration by parts. The interaction in (2.3) can induce a mixing between Ti and 
Z, and as a result the couplings between the Z boson and other SM particles are modified. 
This mixing is constrained by precision electroweak measurements at the permille level 
[61], which implies the bound (with gw the SU(2) gauge coupling) 

cos 9w dh < ]^g-3 

gw 

Note that g]^ is not necessarily suppressed, since My-j can be much larger than mz- Besides 
this mixing effect, the Lagrangian (2.3) does not play a role in our analysis of mono-Z 
production. For the Yq mediator, the leading interaction terms come from 

= Myo 101' lo + gL 101' ^o' • (2.5) 


(2.4) 


When one of the Higgs fields takes its vacuum expectation value, these interactions induce 
a mixing between Yq and the SM Higgs field h, which causes a universal rescaling of all 
Higgs couplings. Current measurements of the Higgs couplings imply the bound [62] 


vMyq ghi 

Myo - K 


< 0.37. 


( 2 . 6 ) 


The interactions (2.5) also give rise to Higgs-boson couplings to DM, which have been 
extensively investigated in the context of Higgs-portal models [63-66] . Since our focus here 
is on the mono-Z signal, we do not take these kind of interactions into account. 

Until now, the interactions between the scalar mediator Yq and SM particles which 
we have considered are either strongly suppressed by Yukawa couplings or do not induce a 
mono-Z signal. This motivates us to include (effective) operators of higher mass dimension. 
The relevant dimension-5 operators are^ 


rYo 

■^EW, dim-5 


A 




!,0\ ' / \ J ' JD fJsi' ' J C. 


(2.7) 


where are the dual field-strength tensors, and A is some large energy 

scale. These effective operators can be induced by loop graphs such as those shown in 
figure 2, where the heavy fermion / is integrated out. The dimension-5 terms in dim -5 
may be considered as a “partial” UV completion of an effective local VVXiXi interaction, 
where the resonance structure of the s-channel Yq mediator can be probed experimentally. 
The underlying physics is that the SM gauge bosons can couple to the mediator Yq through 
loops containing some (very) heavy new particles. Of course, such loop diagrams can also 
induce operators of higher dimension than 5. Their phenomenological impact will however 
be very small provided that the new-physics scale A is sufficiently large. 

^ There exist analogous couplings of Yq to gluons and photons, but they do not contribute to mono-Z 
production at order 1/A. Note, in particular, that the decay chain gg ^ Yq ^ ZZ* —>■ Zvv, which is an 
irreducible background to DM searches in the mono-Z channel, comes with a 1/A^ suppression. 
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Figure 2. Mono-Z production induced by a loop of some heavy fermion /. 



x- 


X 




Figure 3. Different mono-^ production mechanisms at the LHC. Scenarios (a) and (b) are studied 
in our analysis while the scenario (c) is briefly discussed in appendix A. 


The above simplified models are fairly general and flexible. It is straightforward to 
implement them in the universal FeynRules output (UFO) format [67], which allows us to 
study DM production processes within the MadGraph5_aMC@NLO framework in all rele¬ 
vant production channels, within different benchmark scenarios, and with NLO accuracy. 
As mentioned earlier, in the present work we focus on the mono-Z signal for illustration 
purposes, but readers interested in other mono-X signals will find our framework useful as 
a tool that provides accurate and realistic simulations of various DM production channels. 
The UFO files for our simplified models are available at the FeynRules repository [58]. 
Given the list of interaction terms we have specified, the main production mechanisms 
contributing (at tree level) to the mono-Z signal are shown in the first two graphs in fig¬ 
ure 3. The crossed vertices denote the interactions contained in the above Lagrangians, 
while all other vertices are those of the SM. In the type-A scenario shown in the first graph, 
the Z boson is emitted from the initial-state quarks, and the vector mediator Yi connects 
the SM to the dark sector. In the type-R scenario indicated by the second diagram, the 
DM couples to the Z boson through a scalar mediator Yq. In view of the large number 
of relevant couplings that enter the results, we shall in the following discussion choose a 
few representative values for illustration purposes. They are listed in table 1. For the 
masses of DM particles and mediators we follow the suggestions of [40]. A grid with the 
corresponding values is shown in table 2. 

The third diagram in figure 3 shows a topology not covered by our analysis. If the 
sector that couples the SM to the DM consists of very heavy particles, then the resulting 
effective interactions at accessible LHG energies can be represented by contact interactions 
consisting of both visible-sector and dark-sector fields. This gives rise to local, higher- 


- 6 - 




Scenarios 

Dark matter 

Relevant couplings 

Interactions 

Vector mediator Yi 

Ai 

Xd 

= 0.25, = 1 

spin-independent 

A 2 

Xd 

g^ = 0.25, g^^ = 1 

spin-dependent 

A 3 

Xc 

g^ = 0.25, = 1 

spin-independent 

Scalar mediator Yq 

Bi 

Xd 

= 0.25, = 1, A = 3 TeV 

CP-even 

B 2 

Xd 

9w ~ 0.25, gxjy ~ A = 3 TeV 

CP-odd 

B 3 

Xc 

< 7^3 = 0.25, = 1, A = 3 TeV 

CP-even 


Table 1. Benchmark scenarios with representative values of coupling constants in different DM 
production scenarios. In each case, the couplings not shown are set to zero. 


moM [GeV] 

M^ed [GeV] 

1 

10 

20 

50 

100 

200 

300 

500 

1000 

2000 

10000 

10 

10 

15 

50 

100 






10000 

50 

10 


50 

95 

200 

300 




10000 

150 

10 




200 

295 

500 

1000 


10000 

500 

10 






500 

995 

2000 

10000 

1000 

10 







1000 

1995 

10000 


Table 2. DM masses todm and mediator masses Mmed used in the benchmark models. 


dimensional operators, which have been studied, e.g., in [47, 68]. This case is briefly 
discussed in appendix A. In such a scenario our simplified-model framework does not apply, 
because there is no mediator involved. In this case an EFT approach can be useful. Indeed, 
these higher-dimensional operators are still recommended for benchmark studies of mono- 
V signals in [40], because the UV completions of such interactions have not yet been fully 
investigated. 

Before closing this section, we would like to give a comparison between the EFT 
approach and the simplified-model framework for the case of a vector mediator, in order to 
motivate the use of simplified models in our work. In the limit of a large mediator mass, 
i.e. My. ^ 2mDM and My. ^ with q the momentum carried by the mediator, the 

mediator fields can be safely integrated out from the effective theory, resulting in higher¬ 
dimensional operators coupling the SM and DM particles in EFT. However, this limit is 
not always applicable at the LHC. To see how exactly the EFT description approaches the 
simplified models, we compare the total cross sections computed in the two schemes. We 
take the mono-Z production process induced by a vector mediator as a specific example. 
We set the DM mass to 100 GeV and all the relevant couplings equal to 1 for simplicity, 
and we vary the mediator mass between 500 GeV to 10 TeV. In order to investigate the 
dependence of the cross section on the width of the mediator particle, we consider the two 
values Ty-j = My^/3, corresponding to a rather broad resonance, and Ty^ = Myj/(87r), 
corresponding to a narrow resonance [7]. No experimental cuts are applied. The results 
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0.5 1.0 2.0 5.0 10.0 

My, [TeV] 


Figure 4. Ratio of the mono-Z production cross sections at the 8 TeV (red) and 13 TeV (blue) 
LHC, computed in simplified models and EFT. The assumed mediator is a vector Y^. The upper 
and lower limits of the bands correspond to the cases where Ty^ = Myj/(87r) and Ty^ = My^/S, 
respectively. The DM mass is chosen to be moM = 100 GeV. 

are presented in figure 4, which shows that the cross section calculated in simplified models 
is generally much larger than that computed in EFT. The enhancement is due to the s- 
channel resonance, and it is larger in the case of a narrow width. This comparison indicates 
that the EFT does not provide an appropriate framework for predicting DM production 
at the LHC. In order to capture the effects of s-channel resonane, even in cases where the 
resonance mass lies in the few-TeV range, simplified models should be used. 

3 Implementation and validations 

Our computations are performed in the framework of MadGraph5_aMC@NLO [54]. 
The simplified models are implemented in the UFO format [67] by making use of the 
FeynRules package [53]. The one-loop corrections are computed using the MadLoop 
program [69]. It is based on the OPP method [70, 71], in which the ultraviolet divergences 
and rational i ?2 terms are calculated automatically by means of the NLOCT package [72]. 
The infrared subtraction terms for real emissions are generated by MadFKS [73]. The 
matching to the parton shower is performed using the MC@NLO framework [74]. 

We have validated our implementation in several ways. We have calculated the virtual 
QCD corrections for mono-Z production via the scalar mediator analytically and obtained 
results that agree with MadLoop. We have compared the cross sections for the mono¬ 
jet production processes^ induced by vector and axial-vector currents with results from 
MCFM and found good agreement. Finally, by adjusting the couplings and masses of 
specific mediators so as to mimic SM particles, we can compare certain DM production 
processes with corresponding SM processes computed using MadGraph5_aMC@NLO. 
For example, by setting Myi = mz we can compare Yi + Z production with ZZ production 

^Unfortunately, the mono-^ production process has not yet been implemented in MCFM. 
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and Yi + h production with Z + h production, both at LO and NLO. In addition, by setting 
Myo = iTT'h we can compare Yq + Z production with Z + h production. In all cases we have 
found that the results agree within the statistic errors. 

4 Total cross sections 

We start by presenting our predictions for the total cross sections for mono-Z production at 
the 13 TeV LHC in the benchmark scenarios defined in table 1, assuming that the Z boson 
decays to We choose NN23L01 and NN23NLO PDF sets [75] for LO and NLO 

results and use dynamical factorization and renormalization scales, with the default values 
fj,f = = \Ht = \ I where the sum is over all final-state particles. The 

DM and mediator masses are varied over the values shown in table 2. Our parameter choices 
are in agreement with the DM benchmark models suggested in [40] for the early phase of 
the LHC Run-H. The width of the mediator Tmed is determined using the MadWidth 
module [76]. It depends on the masses and coupling constants specific to each benchmark 
scenario, but in general, Tmed/TImed is about a few percent in our parameter choices. In 
this section we do not consider any kinematic cuts. Their effects will be discussed in the 
next section. 

4.1 Benchmark scenarios Aj 

We first consider the benchmark scenarios Ai, A 2 and A 3 , which correspond to simplihed 
models with an s-channel spin-1 mediator Yi. The scenario A 2 is perhaps the most inter¬ 
esting one, since it involves spin-dependent couplings to the DM, for which the LHC can 
set more stringent bounds than direct-detection experiments. The total cross sections and 
iL-factors in benchmark scenario scenario A 2 are shown in table 3 and figure 5, while those 
for scenarios Ai and A 3 are given in appendix B. 

We observe that, for fixed mass of the mediator, the cross section is rather insensitive to 
the DM mass as long as it is below the threshold tudm = Aimed/2, but decreases rapidly for 
heavier DM masses. For a fixed DM mass, the cross section increases when Mmed < 2mDM 
and decreases when Mmed > 2mDM- These are typical features for s-channel resonance 
production, which are not captured in an EFT approach. From appendix B, we also 
observe that the cross sections for benchmark scenario A 3 are approximately one order of 
magnitude smaller than those in scenarios Ai and A 2 . This suggests that it may be more 
promising to search for fermionic DM than for scalar DM. The impact of QCD corrections 
on the cross sections is reflected by the iF-factors, which are defined as the ratios of the 
cross sections computed at NLO and LO. We find that the iF-factors are nearly the same 
for all three scenarios, indicating that the QCD corrections are rather insensitive to the 
coupling structure of the simplified model. The iF-factors decrease as Mmed and tudm 
increase, but mostly stay in a range of roughly 1.3 — 1.5. We thus conclude that NLO QCD 
corrections have a noticeable impact on the mono-Z signal. 

The left plot in figure 5 shows that, for a fixed mediator mass Mmed = 500 GeV, the 
cross section is almost unchanged when moM increases from 1 GeV to 50 GeV, but starts 
to decrease when tudm exceeds 50 GeV. When the DM mass exceeds the threshold Mmed/2 
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todm [GeV] 

[GeV] 

10 

20 

(15) 

50 

100 

(95) 

200 

300 

(295) 

500 

1000 

(995) 

2000 

(1995) 

10000 

1 

8.5 

3.5 

1.0 

0.35 

0.10 

4.5e-2 

1.3e-2 

1.7e-3 

l.le-4 

1.3e-8 

10 

4.6e-2 

5.8e-2 

0.90 

0.34 






1.3e-8 

50 

2.5e-3 


2.9e-3 

6.6e-3 

8.0e-2 

4.1e-2 




1.2e-8 

150 

2.0e-4 




3.0e-4 

8.5e-4 

8.8e-3 

1.6e-3 


l.Oe-8 

500 

3.5e-6 






4.5e-6 

2.8e-5 

7.8e-5 

4.1e-9 

1000 

le-7 







1.4e-7 

1.3e-6 

9.4e-10 

todm [GeV] 

iL-factor 

1 

1.57 

1.46 

1.49 

1.48 

1.42 

1.39 

1.38 

1.35 

1.29 

1.29 

10 

1.49 

1.50 

1.48 

1.47 






1.29 

50 

1.41 


1.42 

1.43 

1.42 

1.41 




1.29 

150 

1.38 




1.38 

1.39 

1.40 

1.36 


1.29 

500 

1.33 






1.34 

1.36 

1.29 

1.23 

1000 

1.21 







1.22 

1.27 

1.09 


Table 3. NLO total cross sections (in pb) and iti-factors for mono-Z production in the channel 
pp Z{-^ at the 13 TeV LHC in benchmark scenario A 2 . We use the short¬ 

hand notation “e-n” for 10“". The values of M„ied shown in parenthesis are used at threshold 

-^med — ■ 



Figure 5. Benchmark scenario A 2 : The pp —>■ 2’(—)■ pApL~) -\-cross section and iL-factor as 
functions of the DM mass (left) and the mediator mass (right). The dashed red line in the right 
plot shows a simple scaling behavior a oc predicted by the EFT approach upon integrating 

out the mediator. The bands represent the scale uncertainties, as estimated by varying the scales 
Pf and pr independently by a factor 2 about their default values. 


for on-shell mediator production the cross section dies out very quickly. It is interesting 
to compare this feature with the results obtained in [42, 43] for the similar case of mono- 


10 

























7 production^ using an EFT approach, where it was found that the cross section drops 
beyond ruDM ~ 100 GeV, but does not become negligibly small until tudm ~ 1000 GeV. 
The reason for the discrepancy is that the invariant mass of the DM pair does not pass 
a threshold if the process is described using an EFT framework. The right plot in the 
figure shows that the cross section decreases as the mediator mass increases. However, 
this decrease is considerably slower than a simple scaling law o oc expected in the 

EFT context. From both plots, we observe an almost constant FT-factor aX K k, 1.4. Note 
that the scale uncertainties are not significantly reduced when comparing the LO and NLO 
predictions, the reason being that the LO prediction is independent of as- 

4.2 Benchmark scenarios Bj 


moM [GeV] 

[GeV] 

10 

20 

(15) 

50 

100 

(95) 

200 

300 

(295) 

500 

1000 

(995) 

2000 

(1995) 

10000 

1 

1.2e-2 

7.2e-3 

2.5e-3 

9.8e-4 

3.3e-4 

1.5e-4 

5.1e-5 

6.9e-6 

3.2e-7 

3.1e-ll 

10 

8.8e-5 

l.le-4 

2.5e-3 

9.7e-4 






3.0e-ll 

50 

6.5e-6 


7.5e-6 

1.6e-5 

3.3e-4 

1.6e-4 




3.2e-ll 

150 

5.8e-7 




8.5e-7 

2.4e-6 

5.1e-5 

6.7e-6 


2.5e-ll 

500 

9.5e-9 






1.3e-8 

8.8e-8 

2.8e-7 

8.6e-12 

1000 

9.7e-9 







l.le-7 

2.9e-7 

8.7e-12 

todm [GeV] 

P-fact or 

1 

1.43 

1.42 

1.39 

1.35 

1.33 

1.30 

1.32 

1.31 

1.17 

1.21 

10 

1.37 

1.37 

1.37 

1.33 






1.17 

50 

1.32 


1.31 

1.32 

1.32 

1.32 




1.29 

150 

1.31 




1.32 

1.32 

1.31 

1.28 


1.22 

500 

1.23 






1.25 

1.27 

1.18 

1.10 

1000 

1.25 







1.26 

1.17 

1.13 


Table 4. NLO total cross sections (in pb) and itT-factors for mono-Z production in the channel 
pp Z{-^ p^pi~) + at the 13 TeV LHC in benchmark scenario Bi. We use the short¬ 
hand notation “e-n” for 10“". The values of M^ed shown in parenthesis are used at threshold 
Mnied = 2 todm- 

We now move to benchmark scenarios Hi, B 2 and B 3 , which correspond to simplihed 
models with an s-channel spin-0 mediator Yq coupling the DM to the electroweak sector of 
the SM. The NLO cross sections and iC-factors for scenario Bi are shown in table 4 and 
figure 6 , while the corresponding results for scenarios B 2 and B 3 are given in appendix B. 
We hrst notice that the cross sections in scenario Bi , and even more so in scenarios B 2 and 
B^, are much smaller than those in the three benchmark scenarios Ai- One of the obvious 
reasons is that the relevant operators are of dimension-5 and thus suppressed by the large 
underlying mass scale, which is chosen to be 3 TeV. One can estimate the total cross sections 
for lower values of the new-physics scale by multiplying them with (3 TeV/A)^, even though 
the actual results will be slightly different as the width of the mediator also depends on 
A. Next, we notice that benchmark scenario B 2 with CP-odd interactions has comparable 

®Prior to our work, there did not exist any NLO results on mono-Z prodnction. 
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Figure 6. Benchmark scenario Bi: The pp — >■ Z{-^ pL^ cross section and iti-factor as 
functions of the DM mass (left) and the mediator mass (right). The dashed red line in the right 
plot shows a simple scaling behavior a oc predicted by the EFT approach upon integrating 

out the mediator. The bands represent the scale uncertainties, as estimated by varying the scales 
fj,f and pr independently by a factor 2 about their default values. 

production cross sections to scenario Bi with CP-even interactions, while scenario P 3 leads 
to much smaller cross sections. This is because the and operators in 

(2.7) each give rise to two powers of momentum in the amplitude, which strongly enhances 
the cross sections, while the {D^{D^(j)) operator in (2.7) leads to a constant factor. 
On the other hand, the P-factors do not show significant differences in the three cases, 
and we again find K ~ 1.3 — 1.5 for most regions of parameter space. 

Figure 6 shows the cross section in benchmark scenario Pi as a function of the DM 
(left) and mediator (right) masses. The behavior is quite similar to that observed for 
scenario A 2 in figure 5, but the drop-off near the threshold seen in the left plot is much 
more steep in the case of Pi. Compared with benchmark scenario A 2 , the PT-factor is 
slightly smaller in the present case, while the reduction of the scale uncertainties at NLO 
is more pronounced. 

5 Kinematic distributions 

The total cross sections presented above display some of the basic features of the mono- 
Z production process at the LHC. However, for practical measurements kinematic cuts 
always play a crucial role. The inclusion of NLO QCD corrections and parton-shower 
effects provide us with a realistic description of differential kinematic distributions. In 
this section, we investigate some important distributions as well as the effects of imposing 
experimental cuts. The parton-shower simulations are performed using PythiaG [77]. 

For each benchmark scenario we choose four different parameter points, (tudMi -^med) = 
(1,10), (1,500), (500,10) and (500,500), all in units of GeV, to present our results. We 


- 12 - 
















10 '’ 




10’^ r 



i LHC 13 TeV 

- NLO+PS, (1,10) 
NLO+PS, (1,500) 

- NLO+PS, (500,10) 
NLO+PS, (500,500) 

- NLO+PS, ZZ 

- NLO+PS, WW 


80 100 

M„ [GeV] 




LHC 13 TeV 

— - LO+PS,(m„^,,M„,^^) 

— NLO+PS, (1,10) 
NLO+PS, (1,500) 

— NLO+PS, (500,10) 
NLO+PS, (500,500) 

— NLO+PS, ZZ 

— NLO+PS, WW 


10-^i 


0 




Q. [GeV] 


Figure 7. Kinematic distributions in the benchmark scenario A 2 for four different parameter 
points. The SM background processes pp —)■ ZZ,WW are also shown. 

have similar distributions as and so we do not show them here. 
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Figure 8. Kinematic distributions in the benchmark scenario Bi for four different parameter 
points. The SM background processes pp ZZ,WW —>■ vupL'^p~ are also shown. pT{pi'^) and 
y{p^) have similar distributions as and y{p,~), so we do not show them here. 
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study the differential cross sections in six observables: the missing transverse energy 
the dilepton invariant mass Mu, the muon transverse momentum and rapidity 

y{p~), and the variables A = and n = —px 'Pt/Pt- observable 

A measures the difference between the missing transverse energy and the transverse mo¬ 
mentum of the Z-boson candidate, which should be small for the associated production 
of DM with a Z boson. The observable D depends on the missing transverse momentum 
and its angle with the direction of the Z-boson candidate. Cuts on these variables are 
very efficient to suppress the backgrounds. For comparison, we show the SM backgrounds 
pp —)• ZZ —>■ and pp —>■ W~^W~ —>■ p,~, which also give rise to a mono-Z 

signal. 

Our results for the two benchmark cases A 2 and Bi are shown in figures 7 and 8, 
respectively. We observe that the IFVF-induced SM background is concentrated in regions 
of small pt{p~) and D and has no resonance peaks in the dilepton invariant mass, and 

hence it should be easy to select proper cuts to suppress this background. On the other 
hand, the ZZ-induced SM background exhibits similar distributions as the DM signal 
processes. From the discussion in previous section, we know that the signals with light 
DM masses and heavy mediators have large cross sections in general because the mediator 
can be on-shell. To search for this kind of signals, stringent cuts on FI™®®, Mu and D 
should be imposed, while other cuts can be kept loose in order to increase the ratio of 
signal over backgrounds. Notice that the signal in the benchmark scenario Bi appears to 
have larger FI™®®, Pt{p~) and D than that in the benchmark scenario ^ 2 - This is because 
the benchmark scenario Bi essentially corresponds to the s-channel production of a Z 
boson in association with the mediator Yq, while the benchmark scenario A 2 essentially 
corresponds to the t-channel production of a Z boson in association with the mediator 
Yi. Because of the denominator in the t-channel propagator, events containing particles of 
high transverse momenta are strongly suppressed. In contrast, the s-channel process has 
no such suppression effect, and thus can have relative larger cross sections in regions of 
high transverse momenta. These three observables can be used as discriminators in case 
any signal is observed. 


6 Discovery potential of the mono-Z signal at the 13 TeV LHC 

In this section we study the discovery potential of the mono-Z signal in the Z —)• pApL~ 
channel at the 13 TeV LHC. In general, the production cross section is a function of the 
relevant couplings, the DM mass and the mediator mass. For illustration purposes, we 
fix the coupling constants according to the suggestions of the ATLAS/CMS DM forum 
and take tudm = 10 GeV, as for light DM the LHC experiments have better sensitivities 
than direct-detection experiments. The mediator mass is varied between 100 GeV and 
1000 GeV. Based on the kinematic distributions studied in the previous section, we first 
impose the basic cuts 

>20GeV, |y(/r±)| < 2.5 (6.1) 
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VVed [GeV] 

Basic cuts 

TT’miss 

Jl/rp 

Mu 

A 

n 

Ccut 

Benchmark scenario A 2 

100 

326 

275 

158 

53.7 

20.4 

0.061 

200 

97.7 

86.6 

59.1 

28.0 

11.7 

0.117 

500 

12.9 

12.0 

9.39 

5.76 

2.73 

0.209 

1000 

1.68 

1.59 

1.32 

0.890 

0.451 

0.265 

ZZ 

4747 

3688 

2101 

1379 

16.0 

3.24 • 10-3 

WW 

988 

479 

82.6 

10.6 

0.487 

4.31 • 10“^ 


Benchmark scenario Bi 


100 

0.966 

0.897 

0.684 

0.466 

0.238 

0.245 

200 

0.331 

0.319 

0.281 

0.228 

0.129 

0.388 

500 

5.09 • 10-2 

5.02 • 10-2 

4.81 • 10-2 

4.35 - 10-2 

2.78-10-2 

0.546 

1000 

6.81 • 10-3 

6.77-10-3 

6.63-10-3 

6.23-10-3 

4.24 - 10-3 

0.622 

ZZ 

4747 

3688 

2101 

1379 

16.0 

3.24 - 10-3 

WW 

988 

479 

82.6 

10.6 

0.487 

4.31 - 10-'^ 


Table 5. Cross sections (in fb) for mono-Z production in the channel pp —>■ Z{—>- 

at the 13 TeV LHC after a series of kinematic cuts, for the benchmark scenarios A 2 (top) and Bi 

(bottom). The last column shows the cut acceptance Ccut- We assume todm = 10 GeV. 


to select the signal events and reduce the SM background. We then apply, step after step, 
a series of more advanced cuts, namely 

> 100 GeV, Mu G [85 GeV, 100 GeV] , A < 0.4, > 80 GeV. (6.2) 

The effective cross sections of the signal and SM backgrounds obtained after applying 
these various cuts are given in table 5 for the benchmark scenarios A 2 and Bi. We can 
see that the backgrounds are efficiently suppressed after the advanced cuts applied. From 
the cut acceptances of the signals, we find that the events with heavier mediator mass can 
pass the selection cuts more easily. More specifically, the cut acceptances in the benchmark 
scenario Bi are higher than those in the benchmark scenario ^ 2 - 

Given the cross sections obtained after applying all cuts, we can estimate the discovery 
or exclusion potentials that can be expected in Run-II of the LHC. The corresponding re¬ 
sults are shown in figure 9. If no signal events will have been detected after accumulating an 
integrated luminosity of 750 fb~^, then the regions with mediator mass Mmed < 1000 GeV 
can be excluded at the 3a level, assuming tudm = 10 GeV in the benchmark scenario H 2 . 
On the other hand, if the mediator mass satisfies Mmed < 890 GeV, then the mono-Z signal 
is going to be discovered at the 5a level before accumulating an integrated luminosity of 
1000 fb~^. The discovery potential of the mono-Z signal in the benchmark scenario Bi is 
also shown in figure 9. Because of the small cross section, the discovery potential is not 
as promising as that for scenario H 2 . In order to show the two cases in one plot, we have 
multiplied the cross section in benchmark scenario Bi by a factor of 9, corresponding to 
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Figure 9. 5a discovery and 3a exclusion limits for DM production at the LHC. If a discovery is 
made, then the regions below the red lines are favored. If no signal is found, then the regions below 
the blue lines are excluded. Results in the benchmark scenario Bi are given for A = 1 TeV. 

lowering the new-physics scale A from 3 TeV to 1 TeV In this case, the LHC with an 
integrated luminosity of 1000 fb“^ can probe or exclude the mediator with a mass lighter 
than about 300 GeV or 400 GeV, respectively, again assuming tudm = 10 GeV. 

7 Conclusions and outlook 

In this work, we have implemented a class of simplified models for DM production via s- 
channel vector or scalar mediators in the FeynRules/MadGraph5_aMC@NLO frame¬ 
work, which allows us to obtain accurate and realistic predictions for DM production rates 
as well as kinematic distributions for Run-II of the LHC. Apart from interactions already 
presented in [55], we have added and validated direct interactions between the mediators 
and the electroweak sector of the SM. Our implementation provides a theoretical basis for 
future DM searches at the LHC. 

As an illustration, we have presented the first NLO QCD predictions for mono-Z 
signals in simplified models including parton-shower effects. We have considered several 
benchmark scenarios including both spin-0 and spin-1 mediators. Our calculation has been 
fully automated and can readily be reproduced. We have presented predictions for the 
total pp ^ Z + production cross sections for the parameter sets suggested by the 
ATLAS/CMS DM forum in [40]. The AT-factors vary in the range of about 1.3 — 1.5, with 
the precise values depending on the DM and mediator masses. This shows that the NLO 
corrections have a noticeable impact on the mono-Z signal and should not be ignored. The 
theoretical predictions of the cross sections become more reliable at NLO and in many 
cases the scale uncertainties are reduced. 

In a next step we have studied various kinematic distributions in order to better un¬ 
derstand the feature of the mono-Z signal. Different coupling structures, DM masses, 

^In the case of Mmed < 500 GeV, the change of A from 3 TeV to 1 TeV has negligible effect on the total 
decay width of Yq since the dominant decay channel is R) —>■ XdXd- 
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and mediator masses can all affect these distributions. From these investigations we have 
obtained important information on suitable cuts needed to select the signal events and 
suppress backgrounds in the search for the mono-Z signal. If a DM signal is observed in 
future, some of these variables may be used to determine the structure of the DM couplings. 
We have also estimated the discovery potential of the mono-.Z^ signal at the 13 TeV LHC 
in several benchmark scenarios. 

The present work can be extended in several ways. First, one should investigate more 
signal channels in the context of simplified models, such as mono-IF, mono -7 and mono- 
h. Different channels are sensitive to different operators in the simplihed models. In a 
combined study one would thus be able to obtain comprehensive understanding of the DM 
interactions. Second, the comparison with present experimental results at the 8 TeV LHC, 
including the mono-V searches as well as dijet resonance searches, would help to constrain 
the parameter space. Third, from our studies we find that the kinematic distributions 
in different scenarios can be rather different. Thus, optimized cuts should be chosen in 
searching for the signals or imposing constraints for specific operators. These interesting 
problems are left to future work. 
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A Contact interactions in EFT 


In the main body of this paper we have focused on simplified models for DM interactions 
with the SM. However, as mentioned in section 2, there is also a case in which the simplified- 
model framework does not apply and an EFT approach should be used. In this appendix 
we provide some results for the particular case where effective local operators describe 
contact interactions between DM particles and the electroweak sector of the SM without 
mediators. Figure 10 shows an example of a loop diagram containing a new, very heavy 
particles / and F', which can generate such operators. They could produce signals with 
electroweak gauge bosons and Higgs bosons in the hnal states [47, 50, 51], and therefore 
should be investigated as well [40]. 



Figure 10. Mono-Z production through a loop-induced operator. / is some very heavy fermion, 
while y is a very heavy boson. 

The complete set of operators up to dimension seven can be found in [47, 51, 68 ]. A 
comprehensive study including the full set of operators is beyond the scope of this work. 
Instead, we only focus on two of them. Our main purpose is to present sample results for 
such loop-induced operators and demonstrate that the approach we have used for simplihed 
models can easily be adapted to the case of EFT studies. In particular, we consider the 
operators 

/:eft = ^ XnXn + W^.W^’^XdXd . (A.l) 

After electroweak symmetry breaking, these operators contain the dimension-5 operator 
XylF+^IT” and the dimension-7 operator both of which are discussed as 

benchmark models in [40]. In fact, the interactions in (A.l) are equivalent to those in 
the benchmark scenarios Hi and H 3 , apart from not having a mediator. Other operators, 
including the dipole operator X^, will not be considered here. 

We have implemented these two operators following the procedure described in sec¬ 
tion 3. We shall consider two benchmark scenarios Ci and ( 72 , corresponding to (cf = 1, 
= 0) and (cf = 0, = 1), respectively. Following the recommendation of [40], we take 

A = 3 TeV and scan over seven DM mass points: = 1,10,50,100,200,400,800,1300, 

all in units of GeV. The resulting cross sections scale like cr ~ A“®, so it is easy to trans¬ 
late our results to lower values of the new-physics scale A. The total cross sections and 
iX-factors are presented in table 6 . We see that for A = 3 TeV the predicted cross sections 
are extremely small and would be impossible to measure at the LHC. However, although 
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ruDM [GeV] 

1 

10 

50 

100 

200 

400 

800 

1300 




Benchmark 

scenario 

Cl 




o-(NLO) [pb] 

1.3e-ll 

1.2e-ll 

8.7e-12 

i.8e-12 

2.9e-12 

8.4e-13 

9.9e-14 

7.8e-15 

K 

1.32 

1.33 

1.31 

1.30 

1.30 

1.24 

1.10 

0.97 




Benchmark 

scenario 

C 2 




(t(NLO) [pb] 

6.2e-8 

6.2e-8 

6.0e-8 

5.5e-8 

4.3e-8 

2.3e-8 

5.4e-9 

7.3e-10 

K 

1.22 

1.23 

1.24 

1.22 

1.21 

1.19 

1.06 

0.91 


Table 6. The NLO total cross sections (in pb) and if-factors for mono-Z production in the channel 
pp —>■ p~) +^T at the 13 TeV LHC in the benchmark scenarios Ci (top) and C 2 (bottom). 

We use the short-hand notation “e-n” for 10“". 



Figure 11. Benchmark scenario C 2 with A = 3 TeV: The pp —>■ Z{-^ p~)+^T cross section and 

AT-factor as functions of the DM mass. The bands represent the scale uncertainties, as estimated 
by varying the scales p / and pr independently by a factor 2 about their default values. 


the energy scale A ~ 3 TeV may thus not be explored, it is still possible to probe lower 
values of A. The difference in the Lorentz structures for scenarios Ci and C 2 leads to sig¬ 
nificant differences in the production cross sections. We also find that the A'-factors in the 
benchmark scenario Ci tend to be somewhat smaller than those in benchmark scenarios 
Ai and Bi. 

In figure 11, we show the cross section in the benchmark scenario C 2 , which is not 
quite as suppressed as the cross section in scenario Ci, as a function of the DM mass. It 
is interesting to compare this figure with the left plot of figure 6. As explained earlier, the 
cross section in EFT drops much slower than that in simplified models. Note also that the 
scale uncertainties are significantly reduced at NLO. Because of the smallness of the signal 
cross section we refrain from showing detailed kinematic distributions for the scenarios 
Cl,2. 

The cross sections of the signal and backgrounds for benchmark scenario C 2 , obtained 
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basic 

TTimiss 

JiL/rp 

Mil 

A 

n 

Ccut 

C 2 

6.13- 10“® 

6.08- 10-5 

5.92 • 10-5 

5.56 • 10-5 

3.82 • 10-5 

0.622 

zz 

4747 

3688 

2101 

1379 

16.0 

3.24 • 10-5 

ww 

988 

479 

82.6 

10.6 

0.487 

4.31 • 10-“^ 


Table 7. Benchmark scenario C 2 with A = 3 TeV: The cross sections (in fb) for mono-Z production 
in the channel pp —>■ Z(— pT) + at the 13 TeV LHC after a series of cuts. The last column 
shows the cut acceptance Ccut- We assume todm = 10 GeV. 



Figure 12. her discovery and 3cr exclusion limits for DM production at the LHC. If a discovery is 
made, then the regions below the red lines are favored. If no signal is found, then the regions below 
the blue lines are excluded. 

after applying the same set of cuts as described in section 5, are given in table 7. With 
A = 3 TeV the cross sections are orders of magnitude too small for detecting a signal at the 
LHC. The corresponding discovery potential and exclusion limits are shown in hgure 12. 
The energy scale in EFT can be constrained to be larger than 650 GeV if no signal is 
observed at the 13 TeV LHC (with an integrated luminosity of 1000 fb“^). This constraint 
does not improve signihcantly even if more data are accumulated. On the other hand, for 
such a low new-physics scale the application of the EFT framework is highly questionable. 
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B Cross sections and i^-factors in benchmark scenarios Ai 3 and B 2,3 

In tables 8-11 we present our results for the total cross sections and A-factors obtained in 
the benchmark scenarios ^1^3 and i?2,3- 


todm [GeV] 

M^ed [GeV] 

10 

20 

(15) 

50 

100 

(95) 

200 

300 

(295) 

500 

1000 

(995) 

2000 

(1995) 

10000 

1 

8.1 

3.6 

1.0 

0.35 

0.10 

4.4e-2 

1.3e-2 

1.7e-3 

l.le-4 

1.3e-8 

10 

8.8e-2 

0.13 

1.0 

0.34 






1.3e-8 

50 

5.1e-3 


6.4e-3 

2.9e-2 

0.10 

4.4e-2 




1.3e-8 

150 

4.7e-4 




8.2e-4 

6.6e-3 

1.2e-2 

1.7e-3 


1.2e-8 

500 

l.le-5 






1.5e-5 

3.5e-4 

l.Oe-4 

6.9e-9 

1000 

4.1e-7 







6.1e-7 

2.1e-5 

2.4e-9 

Todm [GeV] 

if-factor 

1 

1.52 

1.52 

1.52 

1.47 

1.43 

1.39 

1.38 

1.36 

1.30 

1.30 

10 

1.51 

1.50 

1.51 

1.45 






1.28 

50 

1.43 


1.44 

1.46 

1.42 

1.39 




1.31 

150 

1.38 




1.38 

1.41 

1.38 

1.36 


1.26 

500 

1.35 






1.36 

1.37 

1.31 

1.26 

1000 

1.26 







1.27 

1.31 

1.16 


Table 8. NLO total cross sections (in pb) and iti-factors for mono-Z production in the channel 
pp —)■ at the 13 TeV LHC in benchmark scenario Ai. We use the short¬ 

hand notation “e-n” for lO”". The values of Mj^ed shown in parenthesis are used at threshold 
-Mined = 2 todM- 


Todm [GeV] 

M^ed [GeV] 

10 

20 

(15) 

50 

100 

(95) 

200 

300 

(295) 

500 

1000 

(995) 

2000 

(1995) 

10000 

1 

0.98 

0.44 

0.13 

4.2e-2 

1.2e-2 

5.5e-3 

1.4e-3 

1.9e-4 

1.2e-5 

7.9e-10 

10 

2.9e-3 

3.6e-3 

0.10 

4.1e-2 






7.8e-10 

50 

1.6e-4 


1.9e-4 

4.1e-4 

8.2e-3 

4.7e-3 




7.5e-10 

150 

1.3e-5 




1.9e-5 

5.3e-5 

7.4e-4 

1.6e-4 


6.3e-10 

500 

2.2e-7 






2.8e-7 

1.8e-6 

7.4e-6 

2.6e-10 

1000 

6.2e-9 







8.5e-9 

8.3e-8 

6.0e-ll 

rriDM [GeV] 

TL-factor 

1 

1.54 

1.54 

1.52 

1.48 

1.42 

1.38 

1.39 

1.37 

1.32 

1.32 

10 

1.48 

1.49 

1.52 

1.49 






1.33 

50 

1.41 


1.43 

1.42 

1.43 

1.42 




1.28 

150 

1.38 




1.39 

1.38 

1.37 

1.37 


1.29 

500 

1.35 






1.34 

1.37 

1.32 

1.23 

1000 

1.23 







1.24 

1.29 

1.13 


Table 9. NLO total cross sections (in pb) and iL-factors for mono-Z production in the channel 
pp —)■ 2(—>■ pApL~) + at the 13 TeV LHC in benchmark scenario A^. We use the short¬ 
hand notation “e-n” for lO”". The values of M^^ed shown in parenthesis are used at threshold 
-Mmed — 2?nD]vi- 
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TODM [GeV] 

[GeV] 

10 

20 

(15) 

50 

100 

(95) 

200 

300 

(295) 

500 

1000 

(995) 

2000 

(1995) 

10000 

1 

2.0e-3 

1.7e-3 

9.2e-4 

5.3e-4 

2.4e-4 

1.2e-4 

4.3e-5 

6.2e-6 

3.0e-7 

2.8e-ll 

10 

4.8e-5 

6.2e-5 

9.2e-4 

5.3e-4 






2.8e-ll 

50 

7.1e-6 


8.6e-6 

3.0e-5 

2.3e-4 

1.2e-4 




2.8e-ll 

150 

9.2e-7 




1.5e-6 

l.le-5 

4.3e-5 

6.1e-6 


2.6e-ll 

500 

2.2e-8 






3.1e-8 

l.Oe-6 

2.8e-7 

1.2e-ll 

1000 

2.2e-8 







4.3e-6 

2.8e-7 

l.le-11 

todm [GeV] 

iL-factor 

1 

1.38 

1.39 

1.35 

1.33 

1.32 

1.31 

1.31 

1.30 

1.16 

1.23 

10 

1.34 

1.34 

1.34 

1.33 






1.20 

50 

1.31 


1.31 

1.33 

1.31 

1.31 




1.21 

150 

1.31 




1.31 

1.32 

1.31 

1.27 


1.23 

500 

1.25 






1.25 

1.28 

1.18 

1.12 

1000 

1.24 







1.31 

1.18 

1.00 


Table 10. NLO total cross sections (in pb) and if-factors for mono-Z production in the channel 
pp —t Z{^ + J^T at the 13 TeV LHC in benchmark scenario B 2 . We use the short¬ 

hand notation “e-n” for 10“". The values of M^^ed shown in parenthesis are used at threshold 

.^med — ■ 


Todm [GeV] 

[GeV] 

10 

20 

(15) 

50 

100 

(95) 

200 

300 

(295) 

500 

1000 

(995) 

2000 

(1995) 

10000 

1 

1.2e-5 

2.5e-6 

1.2e-7 

9.1e-9 

5.0e-10 

7.6e-ll 

6.0e-12 

1.5e-13 

3.0e-15 

3.3e-18 

10 

1.5e-6 

2.7e-6 

l.le-5 

8.7e-7 






1.7e-16 

50 

1.4e-8 


1.9e-8 

1.6e-7 

l.Oe-6 

1.7e-7 




7.3e-16 

150 

2.6e-10 




5.6e-10 

8.5e-9 

9.9e-8 

2.4e-9 


8.1e-16 

500 

9.8e-13 






1.5e-12 

9.5e-ll 

1.9e-10 

2.6e-16 

1000 

9.8e-13 







5.6e-10 

1.9e-10 

2.6e-16 

Todm [GeV] 

iti-factor 

1 

1.48 

1.47 

1.41 

1.36 

1.33 

1.34 

1.33 

1.39 

1.37 

1.42 

10 

1.43 

1.43 

1.42 

1.35 






1.40 

50 

1.36 


1.35 

1.32 

1.33 

1.33 




1.36 

150 

1.33 




1.33 

1.32 

1.31 

1.31 


1.30 

500 

1.32 






1.31 

1.33 

1.20 

1.30 

1000 

1.31 







1.32 

1.21 

1.27 


Table 11. NLO total cross sections (in pb) and if-factors for mono-Z production in the channel 
pp —)■ Z(—)• pL^p~) + at the 13 TeV LHC in benchmark scenario B^. We use the short¬ 
hand notation “e-n” for 10“". The values of M^ed shown in parenthesis are used at threshold 
.^med — 27nD]V[- 
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